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1.  Background 


In  the  course  of  supporting  many  programs  by  evaluating  inspection  techniques,  the  NDE  group  has 
needed  an  analytical  tool  to  model  or  simulate  the  propagation  of  ultrasonic  waves  through  multi¬ 
layered  structures.  Launch  vehicles  and  spacecraft  have  many  bonded  layered  stractures,  such  as  the 
propellant/insulation/case  assembly  in  solid  rocket  motors,  multilayered  printed  circuit  boards  or  solar 
array/composite  panels  in  satellites.  NonDestructive  Evaluation  (NDE)  techniques  are  used  to  ensure 
bond  integrity  as  part  of  standard  quality  assurance  programs  and/or  to  restore  confidence  when  an 
unusual  event  occurs  that  may  have  compromised  the  hardware.  Ultrasonic  techniques  are  among  the 
most  common  and  powerful  NDE  techniques  used  to  probe  bond  integrity;  however,  beyond  a  single 
layer,  the  signals  can  become  very  complex  as  the  number  of  echoes  grows  with  the  sound  waves 
reverberating  in  each  layer,  and  the  reverberations  spawning  additional  echoes  in  other  layers.  This 
signal  complexity  often  makes  it  difficult  to  determine  inspection  feasibility  without  performing 
expensive  and  time  consuming  tests  on  representative  hardware  samples  with  simulated  flaws.  A 
computer  model/simulation,  however,  provides  initial  feasibility  evaluation  at  a  low  cost  and  in  a 
timely  manner.  Additionally,  the  model  is  a  valuable  tool  for  interpreting  complex  ultrasonic  signals 
in  the  laboratory  and  in  the  field.  The  “ID  Multilayer  UT  Model”  was  developed  under  the  NDE 
MOIE  program  to  address  these  needs. 

“ID  Multilayer  UT  Model”  was  written  in  LabView™  and  is  a  powerful,  flexible,  and  user-friendly 
program.  The  program  allows  one  to  simulate  the  propagation  of  ultrasonic  waves  through  multi¬ 
layered  structures  and  provides  user  feedback  through  real-time  display  of  waves  in  the  layers  and  the 
oscilloscope  trace  that  one  would  expect  to  observe.  In  this  report,  the  Physical  Model  section 
describes  the  physical  principles  behind  the  program;  the  Architecture  section  discusses  the  logic 
employed;  the  Verification  section  compares  simulation  results  with  laboratory  experiment;  and  the 
Using  the  Program  section  describes  the  user  interface.  Detailed  documentation  for  the  LabView™ 
virtual  instruments  is  provided  in  the  appendix. 
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2.  Physical  Model 


A  layered  structure  can  be  inspected  by  using  a  transducer  to  send  sound  waves  through  a  coupling 
medium  in  contact  with  the  top  layer;  see  left  side  of  Figure  1.  As  the  sound  travels  through  the 
structure,  portions  of  the  sound  energy  are  transmitted  or  reflected  at  each  interface.  In  a  pulse-echo 
technique,  the  sound  reflected  back  to  the  top  surface  is  detected  by  the  original  sending  transducer. 

In  the  alternative  through-transmission  technique,  a  second  transducer  on  the  opposite  side  detects  the 
sound  transmitted  through  the  entire  structure.  In  either  case  the  system  in  general  can  be  treated 
analytically  as  n  layers  sandwiched  between  top  and  bottom  semi-infinite  layers.  The  simplifications 
assumed  in  this  model  are  the  following: 


1 .  The  layers  are  flat  and  parallel. 

2.  The  waves  travel  normal  to  the  interfaces.  At  an  interface  a  wave  will  be  reflected  and 
transmitted,  but  not  mode  converted. 

3.  The  material  in  each  layer  is  uniform;  i.e.,  the  ultrasonic  wave  does  not  interact  with 
microstructure  in  a  material  other  than  what  is  accounted  for  by  the  general  attenuation. 

4.  Dispersion;  i.e.,  the  different  frequency  components  of  the  wave  propagating  at  different 
velocities  and  attenuations,  is  not  included. 

5.  Beam  spread  is  not  included,  although,  to  some  degree,  it  can  be  factored  into  the 
attenuation. 


Transducer 


■00§00§0i:.. 


^000: 
0000: 


Physical  Sketch 


mm^m 


Analytical  Model 


Figure  1.  The  physical  sketch  on  the  left  shows  a  typical  ultrasonic  inspec¬ 
tion  of  a  three-layer  structure  with  water  as  the  coupling  medium. 
The  analytical  model,  on  the  right,  is  generalized  to  n  layers  sand¬ 
wiched  between  top  and  bottom  semi-infinite  layers. 
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As  the  ultrasonic  wave  propagates,  two  methods  will  be  considered  to  modify  the  wave — (1)  the  gen¬ 
eral  attenuation  of  the  isotropic  uniform  material,  and  (2)  the  interaction  at  an  interface.  The  wave’s 
pressure  amplitude  is  attenuated  as  the  wave  travels  through  a  material  a  distance  d  by  a  factor  of  e'“ 

,  where  a  is  the  material’s  attenuation  coefficient  in  neepers/mm  or  Np/mm.  Note,  the  attenuation  is 
also  commonly  given  in  units  of  dB/mm,  and  1  Np/m  =  8.686  dB/m.  Typical  attenuation  factors  are 
~5  dB/m  for  A1  and  2.2  dB/mm  for  water  at  1  MHz  and  25°C. 

The  interaction  at  an  interface  results  in  the  wave  energy  partitioning  into  a  reflected  wave  and  a 
transmitted  wave;  see  Figure  2.  The  amplitudes  of  these  waves  are  determined  by  the  reflection  coef¬ 
ficient,  R„  j,^i,  and  transmission  coefficient,  T.  These  coefficients  depend  on  the  two  acoustic 
impedances  of  the  layers,  and  Zn+l-  For  the  interface  between  layer  n  and  layer  n+1,  the  coeffi¬ 
cients  are  the  following: 


Zn+1  + 


r= 


Zw+l  + 


The  waves  are  separated  into  two  arrays;  i.e.,  the  forward  traveling  wave  and  the  reverse  traveling 
wave.  In  accordance  with  the  superposition  principle,  the  total  sound  wave  in  a  layer  is  just  the  sum 
of  the  forward  and  reverse  waves  or  arrays  in  this  case. 


Figure  2.  The  interaction  at  an  interface  between  the  n  and  n+1  layer  with 
acoustic  impedances,  respectively,  of  Z„  and  Zn+l- 
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3.  Architecture 


3.1  Overview 

To  ran  a  simulation,  the  user  enters  acoustic  and  geometric  parameters,  and  required  resolutions.  The 
program  then  proceeds  to  do  the  following:  (1)  calculate  the  array  dimensions;  (2)  request  input 
waveform  information  (at  users  option);  (3)  propagate  the  waves  until  the  user  presses  the  stop  button 
or  pause  button;  and  (4)  display  or  save  to  a  file  the  wave  in  any  layer  and/or  the  scope  outputs 

3.2  Inputs  and  Outputs 

The  inputs  and  outputs  are  shown  schematically  in  Figure  3.  The  user  selects  the  number  of  layers 
between  semi-infinite  top  and  bottom  layers  (Figure  1).  The  required  input  parameters  for  each  layer 
are  the  acoustic  impedance,  Z„;  the  velocity,  v„;  the  path  length  in  layer,  L^;  and  the  attenuation,  (X„. 
Material  parameter  inputs  also  include  the  impedance  of  the  medium  from  which  the  wave  is  input 
into  the  first  layer,  Zj^p,  and  the  impedance  of  the  medium  in  contact  with  the  last  layer,  7^.  The 
user  selects  the  input  wave’s  frequency  and  shape  parameters,  or,  alternatively,  manually  inputs  an 
arbitrary  wave.  The  outputs  include  displays  of  the  waves  propagating  in  individual  layers  and  the 
signals  that  would  be  received  in  pulse-echo  or  through-transmission  techniques.  The  program’s 
received  signal  displays  are  called  Top  Scope  and  Bottom  Scope  since  in  inspection  the  signals  are 
displayed  on  an  oscilloscope  or  scope.  The  Top  Scope  is  the  pulse-echo  signal,  and  the  Bottom  Scope 
is  through-transmission  signal. 


Figure  3.  Schematic  of  program’s  inputs  and  outputs. 
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3.3  Arrays  and  Array  Dimensioning 

The  program  uses  a  number  of  arrays  see  Table  1.  The  arrays  are  either  user  input  values,  internal 
calculation  arrays,  or  output  value  arrays.  As  mentioned  earlier,  sound  waves  are  moving  in  the  for¬ 
ward  and  the  reverse  directions,  and  separate  arrays  are  set  up  in  each  layer  for  each  direction;  see 
Figure  4.  The  dimensioning  of  the  layers’  wave  arrays  is  dependent  on  the  input  wave’s  frequency,  f, 
the  v„’s,  Ln’s,  and  maxAL^’s.  The  input  wave  has  i  points  per  cycle,  and  a  time  between  points.  At 
=l/(f  i).  liie  n*  layer’s  arrays  have  the  number  of  points,  k,,,  calculated  with  the  expression. 


=  Integen 


U„Ar 


y 


However,  to  satisfy  the  user’s  input  resolution  limits,  the  k^’s  and  i  are  calculated  by  an  iterative 
process.  Initially,  i  is  set  to  the  i^^  input  by  user,  and  the  1^’s  are  calculated.  If  the  condition  of  (L„  / 
k„)  <  max  AT  (recall  the  maxAL^’s  are  user  inputs)  is  not  satisfied  for  all  n’s,  then  i  is  increased  until 
the  condition  is  met  for  all  n.  Following  the  calculation  of  k^’s,  the  Prop/el  array  is  calculated.  This 
array  is  the  propagation  attenuation  factor;  i.e.  Prop/el„  =  exp(-a„AL„).  The  TopScope  and  BtmScope 
arrays’  dimensions  are  open-ended  and  grow  with  each  iteration,  i.e..  At,  of  the  model. 


Table  1.  Arrays  used  in  model. 


Array  Type 

Array  Name 

Description 

User  Input 

Z 

velocity 

Atten  dB/mm, 
Length, 

ALmax 

2n,  material  input  impedances 

Vp,  material  velocities 

On,  material  attenuation  coefficients 

Lp,  path  lengths 

maxALp,  maximum  path  differences 

Internal  Calculation 

Input  Wave 

Input  Wave 

(significant  ones) 

FwdWvs 

RevWvs 

pts/Iayer 

RefTop 

RefBtm 

Prop/el 

Forward  Waves  for  all  layers 

Reverse  Waves  for  all  layers 

points  in  each  layer,  kn 

Rp  n-i  at  the  top  interface  of  the  nth  layer 

Rn,n+1  ^t  the  bottom  Interface  of  the  nth  layer 
propagation  values  calculated  from  cx^s 

Out  Arrays 

LayersWvs 

TopScope 

BtmScope 

equal  to  FwdWvs+  RevWvs 
output  signal  from  top  layer 
output  signal  from  bottom  layer 
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3.4  Propagating  the  Waves  (i.e.,  Cycling  Model  One  Iteration) 

One  iteration  of  the  model  corresponds  to  time  being  incremented  by  At,  and  the  waves  in  every  layer 
advancing  one  element.  The  complex  process  of  iterating  the  model  once  in  the  pulse  mode  is  dia¬ 
grammed  in  Figure  4  and  described  by  the  following  array  operations: 

1.  Input  Wave  Array 

th 

a.  Thei  element  is  sent  to  layer  I’s  forward  wave.  All  other  elements  moved  forward 
one  element,  and  a  zero  is  put  in  the  1**  element. 

2.  Layer  1,  Forward  Wave  Array 

a.  The  m*  element  (m  <  kj)  is  multiplied  by  the  attenuation  factor,  Prop/elj,  and  moved 
to  the  m+1  element. 

b.  1  element  is  set  to  sum  of  the  i*  element  from  Input  Wave  and  the  reflected  part 
from  1*'  layer’s  reverse  wave. 

c.  The  original  m  =  kj  element  is  multiplied  by  the  reflection  coefficient  and  sent  to  the 
1®*  layer’s  reverse  wave  and  multiplied  by  transmission  coefficient  and  sent  to  2“ 
layer’s  forward  wave. 

3.  Layer  1,  Reverse  Wave  Array 

a.  The  m*  element  (m  >  1)  is  multiplied  by  the  attenuation  factor,  Prop/elj,  and  moved 
to  the  m-1  element. 

b.  The  kjth  element  is  set  to  sum  of  the  transmitted  part  from  the  2”**  layer’s  reverse 
wave  and  the  reflected  part  from  1**  layer’s  forward  wave. 

c.  The  m  =  1  element  is  (1)  multiplied  by  the  reflection  coefficient  and  sent  to  the  1** 
layer’s  forward  wave  and  (2)  multiplied  by  transmission  coefficient  and  added  to 
TopScope  array. 

4.  Layer  j  with  1<  j  <  n.  Forward  Wave  Array 

a.  The  m*  element  (m  <  kj)  is  multiplied  by  the  attenuation  factor,  Prop/elj,  and  moved 
to  the  m-f-1  element. 

b.  1  element  is  set  to  sum  of  the  transmitted  part  from  the  j- 1  layer  and  the  reflected 
part  from  this  layer’s  reverse  wave. 
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calc 


^^transmitted  part 


Bottom  Scope  Array 


Figure  4.  On  the  left  is  a  general  schematic  showing  the  wave  propagation 
and  interactions  in  the  model.  On  the  right  are  the  details  of  the 
processes  involved  in  iterating  the  model  once  and  moving  the 
waves  (i.e.,  arrays)  one  At. 
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c.  The  original  m  =  kj  element  is  (1)  multiplied  by  the  reflection  coefficient  and  sent  to 
this  layer’s  reverse  wave  and  (2)  multiplied  by  transmission  coefficient  and  sent  to 
layer  j+l’s  forward  wave. 

5.  Layer  j  with  1<  j  <  n.  Reverse  Wave  Array 

a.  The  m*  element  (m  >  1)  is  multiplied  by  the  attenuation  factor,  Prop/elj,  and  moved 
to  the  m-1  element. 

b.  The  kjth  element  is  set  to  sum  of  the  transmitted  part  from  the  j+1  layer’s  reverse 
wave  and  the  reflected  part  from  this  layer’s  forward  wave. 

c.  The  m  =  1  element  is  (1)  multiplied  by  the  reflection  coefficient  and  sent  to  this 
layer’s  forward  wave  and  (2)  multiplied  by  transmission  coefficient  and  sent  to  the 
j-1  reverse  wave. 

6.  Layer  n  (last  Forward  Wave  Array 

a.  The  m  element  (m  <  k„)  is  multiplied  by  the  attenuation  factor,  Prop/el„,  and  moved 
to  the  m+1  element. 

b.  1  element  is  set  to  sum  of  the  transmitted  part  from  the  n- 1  layer  and  the  reflected 
part  from  this  layer’s  reverse  wave. 

c.  The  original  m  =  element  is  (1)  multiplied  by  the  reflection  coefficient  and  sent  to 
this  layer’s  reverse  wave  and  (2)  multiplied  by  transmission  coefficient  and  added  to 
the  BtmScope  array. 

7.  Layer  n  (last  la^er),  Reverse  Wave  Array 

a.  The  m^  element  (m  >  1)  is  multiplied  by  the  attenuation  factor,  Prop/el„,  and  moved 
to  the  m-1  element. 

b.  The  k^th  element  is  set  to  the  reflected  part  from  this  layer’s  forward  wave. 

c.  The  m  =  1  element  is  (1)  multiplied  by  the  reflection  coefficient  and  sent  to  this 
layer’s  forward  wave  and  (2)  multiplied  by  transmission  coefficient  and  sent  to  the 
n-1  reverse  wave. 
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4.  Verification 


The  program  was  experimentally  verified  with  a  simple  three-layer  system  of  polystyrene,  water,  and 
aluminum  in  a  water  immersion  tank.  Figure  5.  A  Panametrics  5052PRX  pulser/receiver  (SN:1482) 
was  used  in  pulse-echo  mode  with  a  Panametrics  20-MHz,  0.25-in.  transducer  (V317/SN:  189384). 
The  waveforms  were  recorded  with  a  Tektronix  720  digital  oscilloscope  (SN:B01 1318)  at  a  digitiza¬ 
tion  rate  of  500  Msamples/s  and  transferred  to  an  Apple  Macintosh  Quadra  950  for  graphing. 

The  program  was  run  using  the  physical  parameters  in  the  table  in  Figure  5  and  with  an  input  wave¬ 
form  of  17.5  Mhz  (the  transducer’s  measured  frequency),  6  half  cycles,  and  a  Gaussian  shape  factor 
of  4.5.  The  overall  comparison  in  Figure  6  shows  very  good  agreement  of  the  simulation  with 
experiment.  The  echoes  are  all  at  the  right  times,  and  the  amplitudes  in  the  early  times  are  consistent. 
The  later  time  amplitudes  of  the  experiment  are  smaller,  which  is  consistent  with  beam  spread  that  is 
not  currently  accounted  for  in  the  program.  The  simulation  shows  a  very  small  echo  at  about  24.2  |j.s, 
which  isn’t  apparent  in  the  experiment  in  Figure  6.  However,  the  echo  is  observable  in  the  experi¬ 
mental  data  at  higher  magnification.  The  difference  in  amplitude  is  probably  due  to  the  standard 
“book”  values  of  impedance  and  attenuation  used  in  the  simulation  being  slightly  off  from  the  prop¬ 
erties  of  the  experiment’s  materials.  The  magnified  portion  of  the  waves  in  Figure  7  illustrates  the 
program’s  accurate  phase  reproduction.  The  phase  flips  on  the  center  echo  in  both  model  and 
experiment. 


Transducer 


■ 


■I* 


Layer 

Material 

L 

mm 

V 

mm/|i,s 

Z 

a 

dB/mm 

@20MH2 

Top 

water 

1.48 

Layer  1 

polystyrene 

14.0 

2.33 

2.46 

0.2 

Layer  2 

water 

3.15 

1.47 

1.48 

0.088 

Layer  3 

aluminum 

9.77 

6.36 

17.46 

0.1 

Bottom 

water 

1.48 

Figure  5.  The  verification  experimental  arrangement.  The  velocity,  imped¬ 
ance,  and  attenuation  values  are  from  the  literature. 
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Figure  6.  The  top  waveform  is  the  Topscope  trace  generated  by  the  ultra¬ 
sonic  model.  It  is  in  very  good  agreement  with  the  scope  trace 
from  the  pulse-echo  experiment  below.  In  the  experiment,  as  the 
wave  travels  farther,  beam  spread  becomes  a  significant  factor  in 
attenuating  the  signal  and  is  not  accounted  for  in  the  model. 
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Figure  7. 
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Ills 

Magnified  region  of  Fig.  6.  The  phase  of  the  model’s  simulated 
signal  matches  that  of  the  experiment  as  seen  in  the  phase  flip  of 
the  reflected  signal  at  ~20.8|is. 


5.  Using  the  Model 


5.1  Overview 

The  major  goals  in  the  interface  design  were  (1)  ease  of  use,  (2)  maximum  flexibility,  and 
(3)  efficient  use  of  screen  space  for  use  on  small  laptops.  The  virtual  instrament  (vi)  panels  are  laid 
out  in  a  compact  logical  format  with  liberal  use  of  color  cues.  The  controls  and  inputs  have  detailed 
description  information  accessible  in  the  LabView™  Help  window.  The  user  has  full  real-time  con¬ 
trol  of  what  results  are  displayed  and  saved.  The  small  and  separate  windows  for  graphs  and  controls 
make  the  program  laptop  friendly. 

5.2  Inputting  Parameters 

The  ID  Multilayer  UT  Model  vi’s  control  panel.  Figure  8,  is  the  primary  method  of  entering  infor¬ 
mation  and  accessing  the  displays  and  secondary  controls.  Many  of  the  values  have  associated  physi¬ 
cal  units;  e.g.,  velocity  has  units  mm/ps.  The  user  is  not  restricted  to  the  default  units  in  displaying  or 
entering  data.  The  units  may  be  changed  by  typing  in  a  consistent  unit;  e.g.,  in/s  can  be  used  instead 
of  mm/|is.  LabView™  automatically  converts  and  keeps  track  of  the  units  in  operations. 


Figure  8.  The  ID  Multilayer  UT  Model  vi’s  control  panel  serves  as  the  pro¬ 
gram’s  main  interface  for  control  and  information  input. 
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Prior  to  running  a  simulation,  the  following  information  is  entered:  (1)  the  geometric  and  acoustic 
parameters  for  the  layers;  (2)  the  resolution  parameters — ^ALmax  for  each  layer  and  #pts/cycle  min; 
(3)  the  number  of  layers;  (4)  input  wave  frequency;  (5)  select  pulse  or  continuous  input  wave;  and  (6) 
select  default  input  wave  or  modify  input  wave.  The  last  two  entries  are  selected  by  buttons  on  the 
panel. 

If  the  choice  was  made  to  modify  the  input  wave,  the  program  pops  up  Waveinput  vi’s  panel.  Figure 
9,  after  the  program  starts  running.  This  vi  either  generates  an  input  waveform  using  the  waveshape 
parameters  input  into  the  program  or  accepts  an  arbitrary  wave  input  point  by  point.  The  input  gen¬ 
eration  parameters  include  the  number  of  half  cycles,  a  Gaussian  exponential  shape  factor,  and  posi¬ 
tive/negative  phase.  The  result  of  a  parameter  change  is  inunediately  observed  in  the  waveform 
graph.  Following  generation  of  an  acceptable  input  wave,  pressing  the  “DONE”  button  exits  Wave- 
Input  vi,  and  the  program  begins  to  wave  propagation. 


Figure  9.  The  input  wave’s  shape  is  controlled  using  the  Waveinput  vi’s  panel. 


16 


The  wave  in  a  layer  or  oscilloscope  output  can  be  displayed  at  any  time  during  the  propagation 
operation.  The  display  vi’s  can  be  opened  by  pressing  a  button  on  the  ID  Multilayer  UT  Model  vi’s 
panel  or  from  a  pull-down  menu.  A  layer  wave  is  displayed  by  a  LayerGraph  vi  panel.  Figure  10,  and 
an  oscilloscope  signal  is  displayed  in  a  UTScope  vi  panel.  Figure  11.  In  the  LayerGraph  vi  panel,  the 
layer  number  control  selects  which  layer’s  wave  is  displayed.  The  unit  button  changes  the  distance 
units  between  millimeters  and  inches.  In  the  UTScope  vi  panel,  a  button  is  used  to  select  pulse-echo 
or  though  transmission  signal  for  display.  The  display  panels  have  on-off  buttons  since  turning  off 
the  graph  displays  can  dramatically  speed  the  propagation  operation. 


Figure  10.  The  wave  in  a  layer  is  displayed  using  the  LayerGraph  vi’s  panel. 

The  layer  number  control  selects  which  layer’s  wave  is  displayed. 
The  unit  button  changes  the  distance  units  between  millimeters  and 
inches. 


Figure  11.  The  pulse-echo  or  though-transmission  oscilloscope  signal  is  dis¬ 
played  in  a  UTScope  vi  panel. 
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Pressing  the  “Save  #s  to  String/File”  button  brings  up  the  Save  #s  Strg/File  vi  panel,  Figure  12.  The 
wave  for  any  layer  and  scope  outputs  can  be  saved  to  files.  The  file  format  is  selected  by  menu  and 
includes  Igor  binary  format  in  addition  to  several  ASCII  text  formats.  The  scope  outputs  may  also  be 
displayed  in  a  string  format  to  allow  copy  and  paste  to  another  application  for  analysis  via  clipboard. 

The  program  will  propagate  the  waves  until  the  “Stop”  button  is  pressed.  All  the  waves  are  still  in 
memory  even  after  the  “Stop”  button  is  pressed  and  can  be  saved  by  running  the  Save  #s  Strg/File  vi. 
During  the  program  operation,  the  scope  arrays  each  grow  by  one  point  for  each  iteration,  thus  the 
program  will  eventually  mn  out  of  memory  if  left  ranning  open  ended.  Currently,  the  program  has 
minimal  error  checking  and  has  no  alerts  for  a  potential  out  of  memory  condition. 


Figure  12.  The  waves  are  saved  with  the  Save  #s  Strg/File  vi. 
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6.  Conclusions  and  Comments 


The  ID  Multilayer  UT  Model  was  verified  with  experiment  to  provide  accurate  simulation  of  oscillo¬ 
scope  outputs.  The  Lab  View™  programming  language  has  eased  the  modification  process  during 
development  and  is  expected  to  ease  the  process  in  future  modifications  and  the  addition  of  new  fea¬ 
tures.  The  program  was  ported  from  Macintosh  to  PC  environment  by  simply  copying  the  vi’s  to  the 
PC  and  opening  them  under  Lab  View™  for  Windows. 

On  the  drawing  board  are  a  few  minor  changes  to  the  program  to  improve  the  ease  of  use.  These 
include  the  following:  (1)  expand  the  Waveinput  wave  generation  options  and  improve  the  arbitrary 
wave  input  interface;  (2)  enhance  the  error  checking;  (3)  expand  the  number  file  formats  for  saved 
wave  data,  and  (4)  add  a  beam-spread  algorithm.  A  more  detailed  manual  is  planned  that  will  docu¬ 
ment  some  of  the  convenient  features  not  in  this  report,  such  as  the  little  database  vi  that  contains 
many  common  materials’  acoustic  properties. 
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Appendix 


ID  Multilayer  UT  Model’s  Virtual  Instruments 
Lab  View™  Documentation 
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THE  MAJOR  VIRTUAL  INSTRUMENTS’  FUNCTIONS 


Page  1 


1D  MultiLayer  Model 

Last  modified  on  10/7/96  at  3:20  PM 

Printed  on  10/7/96  at  3:21  PM 

Connector  Pane 


Q540ppcHD:LabViEW4ppc:UTModel2.llb:1D  MultiLayer  Model 

VI .0,  6  August  1996  by  Jim  Lhota.  Ultrasonic  ID  finite  element  model  for  wave  propagation  through  multilayer 
stack  of  materials.  Created  by  Jim  Lhota  based  on  my  previous  2  layer  model.  The  number  of  layers  is  limited 
currently  to  255,  but  can  easily  be  increased  by  changing  from  U8  integer  to  U16  or  U32  integer. 


Controls  and  Indicators 


Ztop 

Accoustic  Impedance  of  the  semi-infinite  medium  above  the  first  layer  -  (mm/ps)*(g/cm^3)  or 
0.001  *g/(mm^2  ps)  or  (10^6)*kg/(m'^2  s) 

Zbottom 

Accoustic  Impedance  of  the  semi-infinite  medium  below  the  last  layer  -  (mm/ps)*(g/cm^3)  or 
0.001  *g/(mm''2  ps)  or  (10‘^6)*kg/(m''2  s) 

Length 

Thickness  of  each  layer,  i.e.  path  length. 


velocity 

Velocity  of  each  layer. 


Z 

Accoustic  Impedance  of  each  layer  -(mm/ps)*(g/cm^3)  or  0.001  *g/(mm'^2  ps)  or  (10^6)*kg/(m''2 


Atten  dB/mm 

Attenuation  coefficient  of  each  layer. 


25 


1D  MultiLayer  Model 

Last  modified  on  10/7/96  at  3:20  PM 

Printed  on  10/7/96  at  3:21  PM 
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i  [SGL] 


[  [SSL]  I 


ALmax 

Limits  of  maximum  size  of  an  element  in  each  layer.  ALmax  >  Length/#pts. 


#pts/cycle  min 
Minimum  pts  in  a  wavelength. 

CW/Pulse 

Pulse  -  the  waveform  is  input  followed  by  zeros. 

Continuos  -  the  waveform  is  repetively  input  (note:  for  continues  use  an  even  number  of  cycles) 
Frequency 

Wvinput  Ctl 

Brings  up  "WavelnpuLvi"  control  panel  to  set  wave  parameters. 

Graph  Update  Iteration 

stop 

Stops  the  calculation  and  exits.  Note  after  stopping  the  last  data  is  still  in  the  global.vi  and  can  be  saved 
by  running  the  subVI  "Save  #s  to  Strg/File" 

Pause 

fLayers 

The  number  of  layers  to  be  used  in  calculation.  Other  layers  will  be  ignored. 

Save  Waves 

Scopel 

Scope2 

Graphi 

Graph2 

Pts/Layer 

Number  of  elements  that  each  layer  is  divided  into  for  calculation. 

Iteration 

RefBtm 

Calculated  reflection  coefficient  at  the  bottom  of  layer  n  =  (Zn+1  -  Zn)/(Zn+1  +  Zn).  Note  for  last  layer 
Zn+1  =  Zbottom 

RefTop 

Calculated  reflection  coefficient  at  the  bottom  of  layer  n  =  (Zn  -  Zn-1)/(Zn  +  Zn-1).  Note  for  first  layer 
Zn-1  =  Ztop. 

Prop/el 

Propagation  coefficent  between  elements  in  each  layer.  Calculated  from  attenuation  coefficient.  Note 
this  value  is  <  1 ,  e.g.  if  no  attenuation  the  whole  signal  is  propagated  and  the  value  would  be  1 . 

Irssi  il 
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1D  MultiLayer  Model 
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[ui6]  Arrindx 

-  The  waves  for  ali  iayers  are  in  two  giant  arrays  (WvsFwd  &  WvsRev)  and  this  index  is  where  each 

iayers  waves  start  in  these  arrays. 


j  [sail  I  Prop  n_n+1 

^  *  Propagation  coefficent  between  elements  between  layer  n  and  n+1  based  on  the  attenuation  coefficient 

and  not  to  be  confused  with  the  transmission  coefficient  from  the  impedance.  Note  this  vaiue  is  <  1 ,  e.g. 
if  no  attenuation  the  whole  signal  is  propagated  and  the  value  would  be  1 . 
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1D  Multilayer  Model 

Last  modified  on  10/7/96  at  3:20  PM 
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1D  MultiLayer  Model 
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Position  in  Hierarchy 
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List  of  SubVIs 
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ID  MultiLayer  Model 
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Find 

First 

Error 


1D  MultiLayer  Model 

Last  modified  on  10/7/96  at  3:20  PM 

Printed  on  10/7/96  at  3:21  PM 
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Calc  Arrays'  Sizes 

Q540ppcHD:LabVlEW4ppc:UTModei2.llb:Caic  Arrays'  Sizes 
llayer  move  wvs  leiem.vi 

Q540ppcHD:LabVIEW4ppc:UTModei2.ilb:1  layer  move  wvs  leiem.vi 
Interface  Reconcile. vi 

Q540ppcHD:LabVIEW4ppc:UTModei2.llb:lnterface  Reconcile.vi 

Global. vi 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:Globai.vi 

UTLayerGraphI 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:UTLayerGraph1 

UTLayerGraph2 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:UTLayerGraph2 

UTScopel 

Q540ppcHD:LabVIEW4ppc:UTModel2.ilb:UTScope1 

UTScope2 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:UTScope2 
Save  #s  to  Strg/Flle 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:Save  #s  to  Strg/File 
Acoustic  properties 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:Acoustic  properties 
Waveinput.vi 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:Wavelnput.vi 
Close  Panel.vi 

Q540ppcHD:LabVIEW4ppc:vi.lib:Utility:victi.llb:Close  Panei.vi 

Open  Panel.vi 

Q540ppcHD:LabVIEW4ppc:vi.lib:Utility:victl.ilb:Open  Panei.vi 


History 


“ID  MultiLayer  Model  History" 
Current  Revision:  35 
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Calc  Arrays'  Sizes 

Last  modified  on  10/7/96  at  3:59  PM 
Printed  on  10/7/96  at  4:00  PM 
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Connector  Pane 


No.  Layers 
Frequency 
#pts/cycle  (minimum) 


Length 
ALmax 
velocity 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:Calc 


#pts/cycle 

ArrayDim 


Arrays'  Sizes 


Calculates  array  sizes  for  UT  1-D  model  of  multiple  layers  of  materials  sandwiched  between  two  semi-infinite  media. 
Jim  Lhota  July  1996 

Front  Panel 


#pts/cycle 

Frequency 

Length 


velocity 


(minimum) 


ALmax 


[U16] 


No.  Layers 

#pts/cycle 

ArrayDim 
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Calc  Arrays'  Sizes 

Last  modified  on  10/7/96  at  3:59  PM 
Printed  on  10/7/96  at  4:00  PM 


Block  Diagram 


Position  in  Hierarchy 


List  of  SubVIs 


History 


"Calc  Arrays'  Sizes  History" 
Current  Revision:  50 
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Wavelnput.vi 

Last  modified  on  10/7/96  at  4:00  PM 
Printed  on  10/7/96  at  4:00  PM 

Connector  Pane 


#Half  Cycles 
#pts/cycle 
Exponential  Factor 
Done 

Q540ppcHD:LabVIEW4ppc:UTModel2. lib:  Wavelnput.vi 


Front  Panel 


Controls  and  Indicators 

#pts/cycle 
Gen/Input  Wv 
User  Input  Wv 

|f^fl 

#Half  Cycles 
Exponential  Factor 
phase 
Done 


Wave 

isifij 

Last  pt 


inntii 

|gPI 


TF 


im 


I  [SSL] 


U16 
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Waveinput.vi 

Last  modified  on  10/7/96  at  4:00  PM 
Printed  on  10/7/96  at  4:01  PM 
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Block  Diagram 


o  o  n 


Position  in  Hierarchy 


List  of  SubVIs 


Global. vi 

Q540ppcHD;LabVIEW4ppc:UTModel2.llb:Globai.vi 


History 
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Waveinput.vi 

Last  modified  on  10/7/96  at  4:00  PM 
Printed  on  10/7/96  at  4:01  PM 

"Waveinput.vi  History" 

Current  Revision:  40 
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1  layer  move  wvs  lelem.vi 

Last  modified  on  10/7/96  at  4:01  PM 

Printed  on  10/7/96  at  4:01  PM 


Connector  Pane 


RefITopInterface 
FwdWvJn 
Trans/element 
RevWvJn 

RefIBtmInterface 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:1layer  move  wvs  lelem.vi 

Inputs:  Inputs  from  rt  &lft  interfaces,  Afwd_in(i),  Arev_in(i),R  reflection  coef  for  left  and  right  interfaces, 
T=transmission  per  element  of  travel  (i.e.  1-atten). 

1.  Afwd_out(i)  =  T*AfwdJn(i-1)  i^^O; 

Afwd_out(0)  =  Rlftt*T*Arev_in(0)+lnLftlnterface 

2.  Arev_out(i)  =  T*Arev(i+1)  i^tn; 

Arev_out(n)  =  Rlft*T*Afwd(n)+lnRtlnterface 

3.  OutLeftlnterface=Arev_in(0)-Rlftt*Arev_in(0) 

4.0utRightlnterface=Afwd_in(n)-Rlftt*AfwdJn(n)s 

Front  Panel 


1 

FwdWv_out 
OutTopInterface 
OutBtmInterface 
j  L=~  RevWv_out 


l^yi 


FwdWv_in 


RevWv_in 

Trans/element 

RefITopInterface 

RefIBtmInterface 

FwdWv_out 

RevWv_out 
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llayer  move  wvs  lelem.vi 

Last  modified  on  10/7/96  at  4:01  PM 

Printed  on  10/7/96  at  4:01  PM 

fiinj 

[i’iSjl]  OutTopInterface 
|j  SCL  jl  OutBtmInterface 


Biock  Diagram 


List  of  SubVIs 


History 

"llayer  move  virvs  lelem.vi  History" 
Current  Revision:  48 
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interface  Reconcile. vi 

.ast  modified  on  10/7/96  at  4:02  PM 

Printed  on  10/7/96  at  4:02  PM 

Donnector  Pane 


n_0  ut  Btm  I  nte  rf  ace 
n_RevWv_in 
PropCoef_n_n+1 /element 
n+1_FwdWvJn 


n_RevWv_out 

n+1_FwdWv_out 


n+1_OutToplnterface - 

Q540ppcHD:LabVIEW4ppc:UTModel2.llb:lnterface  Reconcile. vi 


Inputs:  Inputs  from  rt  &lft  interfaces,  Afwdjn(i),  ArevJn(i),R  reflection  coef  for  left  and  right  interfaces, 
T=transmission  per  element  of  travel  (i.e.  1-atten). 

1.  Afwd_out(i)  =  T*Afwd_in(i-1)  i^^O; 

Afwd_out(0)  =  Rlftt*T*ArevJn(0)+lnLftlnterface 

2.  Arev_out(i)  =  T*Arev{i+1)  i?in: 

Arev_out(n)  =  Rlft*T*Afwd(n)+lnRtlnterface 

3.  OutLeftlnterface=Arev_in(0)-Rlftt*Arev_in(0) 

4.0utRightlnterface=Afwd_in(n)-Rlftt*Afwd_in(n)s 


Interface  Reconcile.vi 

Last  modified  on  10/7/96  at  4:02  PM 

Printed  on  10/7/96  at  4:02  PM 

^rssL  fi 


[ssi]  I  n+1_FwdWv_out 


Block  Diagram 


List  of  SubVIs 


History 

"Interface  Reconcile.vi  History" 
Current  Revision:  58 
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UTScopel 

Last  modified  on  10/7/96  at  4:03  PM 
Printed  on  10/7/96  at  4:03  PM 


Connector  Pane 


Q540ppcHD:LabVIEW4ppc:UTModel2.ilb:UTScope1 

Displays  scope  trace  from  ultrasonic  models.  VI.O  Jim  Lhota,  1  Aug  1996. 
Front  Panel 


Controls  and  Indicators 


|l  TF  il  Graphi 
stop 

j  [&Si]  I 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security  programs,  spe¬ 
cializing  in  advanced  military  space  systems.  The  Corporation's  Technology  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staffs  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  Technology  Centers: 

Electronics  Technology  Center:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  Micro-Electro-Mechanical  Systems  (MEMS),  and  data  storage 
and  display  technologies;  lasers  and  electro-optics,  solid  state  laser  design,  micro-optics, 
optical  communications,  and  fiber  optic  sensors;  atomic  frequency  standards,  applied 
laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam  control, 
LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of  new 
materials:  metals,  alloys,  ceramics,  polymers  and  composites;  development  and  analysis 
of  advanced  materials  processing  and  deposition  techniques;  nondestructive  evaluation, 
component  failure  analysis  and  reliability;  fracture  mechanics  and  stress  corrosion;  analy¬ 
sis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  spacecraft  structural 
mechanics,  space  environment  effects  on  materials,  hardening  and  vulnerability  assess¬ 
ment;  contamination,  thermal  and  structural  control;  lubrication  and  surface  phenomena; 
microengineering  technology  and  microinstrument  development. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing 
using  atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analy¬ 
sis;  effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth's  atmos¬ 
phere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate  radia¬ 
tions  on  space  systems;  space  instrumentation;  propellant  chemistry,  chemical  dynamics, 
environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of  mis¬ 
sile  plumes,  and  sensor  out-of-field-of-view  rejection. 


